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In this letter, the problem of radiation in a fiber geom¬ 
etry interacting with a two level atom is mapped onto the 
anisotropic Kondo model. Thermodynamical and dynamical 
properties are then computed exploiting the integrability of 
this latter system. We compute some correlation functions, 
decay rates and Lamb shifts. In turn this leads to an analysis 
of the classical limit of the anisotropic Kondo model. 

PACS numbers: 42.81.Dp, 72.10.Fk, 75.10.Fk. 


I. INTRODUCTION 

The interaction of radiation with two-level atoms is a 
fundamental problem with many important applications. 
The model goes under the name of the Maxwell-Bloch 
theory in quantum optics. The problem has been ex¬ 
tensively studied theoretically in a framework where the 
radiation is essentially classical and the so-called slowly 
varying envelope and rotating wave approximations are 
commonly made, which are valid near resonance [Q. In 
one dimension, fully quantum integrable versions of this 
system, refered to as the reduced Maxwell-Bloch theory, 
have been solved exactly |^. In some applications, off- 
resonance effects are likely to be significant; for example 
in doped photonic band gap materials. In this letter, we 
show that the single atom case is exactly solvable without 
making any near-resonance approximations, by mapping 
the problem onto the anisotropic Kondo model. Further¬ 
more, we are able to compute important physical quanti¬ 
ties that have hitherto been inaccessible in previous treat¬ 
ments of even the approximate reduced Maxwell-Bloch 
theory. In particular, we compute thermodynamic prop¬ 
erties in the form of impurity corrections to the Stephan- 
Boltzman law, and some electric field correlators. From 
the latter we obtain the atomic spontaneous emission de¬ 
cay rate and the Lamb shifted energy splitting. 


II. THE MODELS 


The interaction of radiation with a two-level atomic 
impurity at Xg is described by the total hamiltonian 


^ _ ^atom _j_ ^field _|_ jjint 


( 1 ) 


with 


1 Z*®® 

dx[[dtcff + (2) 

jjatom ^ ^ , 

and the conventions [ag, a±] = ±2a±, [cr+, (t_] = 173. We 
have made a reduction from the Maxwell theory to a one¬ 
dimensional fiber geometry, following the conventions in 
1^. The dimensionless coupling constant j3 is determined 
by the strength d of the dipole transition and the effec¬ 
tive cross-sectional area of the fiber, /3 = ■sjd. The 

coupling /3 is also related to the lifetime t of the excited 
state of the atom; to lowest order, 7 = 1/7 = /3^wo/4. 

The hamiltonian is known in the condensed matter 
literature as the spin boson hamiltonian [Q, |^. It can 
be mapped onto the anisotropic Kondo hamiltonian Hk- 
dehning the unitary operator 


C/ = -^ exp ( i^a 3 (j){xo) ) (cts -b cr+ -b cr_) 


then H = WHkU, where 


Hk = 


cr_e 




( 3 ) 

( 4 ) 


This allows us to formulate the quantum mechanics in 
either the “optical picture” based on H or the “Kondo 
picture” based on Hk- 

In order to complete the map to the conventional 
Kondo model on the half-line a: < 0, we fold the system. 
Let 4>L{x-\-t), (j)}i{x — t) be the left and right-moving com¬ 
ponents in the bulk, with = 1/7 -b(/_R. Let us set Xg = 0, 
and define even and odd helds in the region a: < 0 


‘pKx, i) = /3 ((/>l(x, t) + (fni-x, t)) 
t)= (3 {(I)l{x, t) - (fni-x, t)) 

(p^ji{x,t) = ipl{-x,t), ip°K[x,t) = -ipl{-x,t). (5) 

Introduce new fields -b The odd and 

even fields decouple; the odd held just becomes a free 
held, while for the even held the hamiltonian reads 

Hr = l [ dx -b 

^ J—00 ^^9 

+ ^(^^e*0(o)V2 + ^_g-.r(o)/2)^ (g) 

where g = 0 ^/Att. 

Recent progress in the understanding of the anisotropic 
Kondo problem will now allow us to compute exactly the 
quantities of interest in the optics problem. 
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III. THERMODYNAMICS 


( 11 ) 


We start with thermodynamics. Consider the sys¬ 
tem at finite temperature T. The partition function 
of the optical system is identical to the Kondo one: 
Z = Tre~^^^ = . Without the atom, the 

photon energy spectrum is given by the usual Plank dis¬ 
tribution. The total partition function is 


T /■“ 

logZ|g=o = -^ j dAlog(l-e-l'=l/^)+logZa 


= -LT + 

6 


log(2coi,h(^)) 


( 7 ) 


where L is the length of the system. Of particular interest 
is the mean energy density 


Cf(T) = C<^logZ = ^r"-^ta„h(^). (8) 

The first term in (^) is the Id analog of the Stephan- 
Boltzman law. 

In the interacting theory, we work with equation (|^). 
The partition function factorizes into a bulk term and 
an impurity contribution Z = 2'buik • -^imp- The bulk 
contribution, including both is the same as (|^): 

log Zbuik = ^LT. 

The impurity contribution in the context of the Kondo 
problem has recently been studied using novel techniques 
in i0- In particular, the partition function is the 
trace of the monodromy matrix studied in detail in Q. 

It is convenient to express physical quantities in terms 
of a “Kondo temperature” Tk, defined such that the lead¬ 
ing low temperature behaviour is logZimp ^ Tk/T. The 
Kondo temperature is then a fixed function of the energy 
splitting and g 



The leading high and low temperature behaviors for ar¬ 
bitrary g are then as follows. At low temperatures, 

\0gZimpi< ^ +J2<^n{g) ( 10 ) 

n—0 ^ 


where ^ tan ’ 


^For the conventions in Q , the ‘spectral parameter’ argu¬ 
ment A of the monodromy matrix T{X) is a function of uio and 
T, namely A = tJo(27rT)®“^/2. The parameter uo has anoma¬ 
lous scaling dimension —g, which is consistent with the the 
renormalization group equation derived in the optical picture 
in ||. 


Oi 


7r2(l-45)(4-5) i 

f l-2g) 
^2-2g) 

iri 

( 2 - 29 ) 

360 p3 1 

( 2-3g) 
( 2-2g J 

iri 

UAJ 


Higher an{g), up to n = 3, can be deduced from results 
in ||]. At high temperatures, one finds 


Ziuip — 2 ■ 


1 ( 5 ) 


Ie. 

T 


2-2g 


O 


T 


where 


T 


( 


2-3g ^ 

2-2g, 


= 1 1VE!£( 

ryi-g) 


4-4g 


2-2g 


( 12 ) 


(13) 


The mean energy density U{T) has the form U = 
C^buik + C^imp where t/buik = 7rT^/6. At low and high 
temperatures for arbitrary g, one finds 


U; 


imp 




T2\ 

’ 


T /Tj 


-a(l-5)T(^ 


2-2g 


T^Tk, 


T > Tk. (14) 


In the quantum optical context, g is generally very 
small, thus it is more relevant to expand Zi^p in powers 
ofg. For small g, Tk ~ i47o/2. In [^, the low temperature 
expansion of Zi^p is expressed in terms of the so-called 
local integrals of motion En-i- With the conventions of 
1^, one can show that 

9"-'b»-.(9)l,=0 = (15) 


where B 2 n are the Bernoulli numbers. Using this, one can 
show that to order g the partition function at arbitrary 
temperature is given by 


^imp = 2 cosh 


Re 


Ie 

T 

4- ( 1-h 


25 ^ sinh 


iTk \ _ 
ttT ) 


log 


Ie 

T 

Ie 

ttT 


(16) 


where 47 ( 2 :) = 9a;logr(a:). One can easily check that 
( 0 ) correctly reproduces the high and low temperature 
expansions in & & at small g. 

The classical limit of the sine-Gordon model thermo¬ 
dynamics had been partially studied in [^. We see here 
that there are, in fact, physical probes exploring the dif¬ 
ferent Bethe ansatz excitations: these are the spin j im¬ 
purities, whose thermal properties depend on the pseu¬ 
doenergies 62 j of the thermodynamic bethe ansatz. For 
instance, when 5 = 0, Zi^p corresponds to the partition 
function of a decoupled two-level system with energies 
±u}o/ 2 : Zimp = 2cosh(u)o/2T), and this can be recovered 
using the expressions for ei in 

The mean energy density can be easily obtained for all 
T from and (fi)- 


2 



















( 22 ) 


IV. PHOTON CORRELATION FUNCTIONS 


Let E denote the relevant component of the elec¬ 
tric field perpendicular to the fiber. Given some ini¬ 
tial state |i), or some mixed state with density matrix 
p = one is interested in the intensity of the 

radiation in the electric field 

Hx) = (17) 

where {0)p = TrpO, and E^~'> {E^^'>) is the creation 
(annihilation) part of the field, E = E^'^'> -\- E^~'>. In the 
quantum measurement theory of Glauber, I{x) is pro¬ 
portional to the number of photons detected at x given 
the initial state specified by p. One is also interested in 
the power spectrum S{uj) of detected photons satisfying 
duiS{uj) = I{x). One can express the above quanti¬ 
ties in terms of the usual electric field correlators. Taking 
as the definition of E^^'^: 

7 r°° ^ 

= ( 18 ) 

one can show that if \ijj) is an eigenstate of the hamil- 
tonian with energy E.^, then {iIj'\E^^'>\iIj) = {tp'\E\ilj) if 
E^' < Eip and zero otherwise. Likewise, {ijj'\E^~'>\il}) = 
{ip'\E\ijj) if E-,p/ > E-,p and zero otherwise. Inserting com¬ 
plete sets of states in & , one finds 

S{uj) = J dTe~"‘^'^{E{x,T)E{x,0))p. (19) 


The power spectrum S{u!) is non-zero only for den¬ 
sity matrices p containing excited states, e.g. a ther¬ 
mal distribution, since the integral in (H) is proportional 
to 9{—uj) for p = |0)(0|. For example, such correlators 
were computed for the harmonic oscillator version of the 
present problem in [pd| . In this paper, in order to illus¬ 
trate our approach to correlation functions, we deal with 
the simpler quantity (0|iil(a:,t)i?(a:,0)|0), which probes 
how easy it is to add a photon to (rather than detect a 
photon in) the vacuum. In the reduction to one dimen¬ 


sion, E = ^ 


Thus we define 


poo 

i(x) = (0| : dt</>(x,t)dt(p(x,t) : |0) = / duj s(uj). (20) 

Jo 

We define the normal ordering in (|2^) such that the free 
photon vacuum energy |w|/47r is subtracted from s(w). 
To lowest order in perturbation theory in the optical pic¬ 
ture, one has 


s(w) 


^ gUJo UJ _ ^2i\uix\^ 

2 Lu'^ — loq + ie 


( 21 ) 


All the previous quantities were defined in the unro¬ 
tated and unfolded theory. The rotation by the unitary 
operator U changes the operator dt(j) to 


dt4>{x,t) dt(j){x,t) - ^asS^x). 

Since the detector is taken to be away from the impurity 
(and at a: < 0 in our conventions), the second term of the 
previous expression can be dropped. Then under folding, 
we find 

dt(t>{x, i) = {dtctt'-ix, t) + dt(j)°{x, t)) , (23) 

where now the odd and even field are defined on the 
negative axis. The odd field is free and has Dirichlet 
boundary conditions, 4>°{0) = 0 , whereas the even field 
has Neumann boundary conditions and interacts with the 
boundary. When computing s{io) the odd field will not 
contribute because the Wick contraction is substracted. 
Thus 

= (24) 

Ibirg 

In order to compute these correlation functions, we 
will make use of the techniques used for the anisotropic 
Kondo model in |^. There, the “natural” basis of the 
excitations of the system, as described by the massless 
sine-Gordon model, is used to compute correlation func¬ 
tions. At first sight, the problem appears untractable 
since, as 5 ^ 0, the spectrum of the sine-Gordon the¬ 
ory consists of a very large number of particles: soli- 
ton, anti-soliton and n bound states or breathers where 
n ~ l/g. The problem simplifies however when con¬ 
sidering the different contributions of these particles to 
the correlators. A quick way to see this is to consider 
< d^(j){z,z)dz'4'{z\z') >= 2gj{z — z')'^, which is unaf¬ 
fected by the boundary interaction. (Above, z = x + iy, 
1 = x — iy, and y = it) This correlator can also be com¬ 
puted using form factors, and expands as a sum of pos¬ 
itive contributions. This gives a good idea of the order 
of magnitude of neglected terms when only the most im¬ 
portant excitations are taken into account: for instance 
the one breather is exact to order g and is an excellent 
approximation for values of g up to g ^ 1/5 (see below). 

The general expression for the current-current correla¬ 
tor is 

cxD ^ n n 

< dz(i){z,z)dz'(i){z',z') >= XI XI / n ^ (2^) 


where is the form factor of the right moving current 
operator and the dot product means that there is a con¬ 
traction between the isotopic indices of the R matrices 
and that of the form factors. 

The one breather has mass Mi = 2 sin( ^)(^^ ) and using 
the results of we find the intensity to be 


7 ^ r ^ 


47r 


A:i(log( 


MiTi 


-))-l 


^ — 2ujx 


(26) 
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where for convenience we have defined Tb = 
TR:/tan( 7 r( 7 / 2 (l — gj), and Ki is related to the reflection 
matrix for the breather by Ki{6) = Ri{i^ — 9), where 


Ri[e) = 


tanh(f 

tanh(| 


(27) 


The ‘power spectrum’ s(a;) thus takes a very simple form 

, ,2 


■S(^) = 

where 

So = 2 cos 


27r y — Sq + 


■Kg 


2-2g 


Tk, 7 = 2 sin 


1-5 


(28) 

Tk- 

(29) 


This result reproduces the perturbation theory compu¬ 
tation at small g dUl)- Near resonance oj « Sq, s(u;) 
takes on the familiar Wigner-Weisskopf Lorentzian line- 
shape, with a Lamb shifted resonant frequency So and a 
half-width j/2: 


s{uj) 


T f _\ —2iujx 

8k V (w - So)2 -I- ( 7 / 2 )^) 


(30) 


We remark that since 7 is twice the imaginary part of 
the pole of i?i in the momentum variable k = MiTee®, 
in a sense 7 represents the exact decay rate into the first 
breather. In the small g limit, the first breather is iden¬ 
tified with the photon and thus 7 agrees with the lowest 
order perturbative result 7 = KgujQ. 

One can check that this one-breather approximation is 
better than the latter perturbative result. Actually, it 
turns out that it is an excellent approximation for the 
correlators for values of g at least up to g = 1/5! Con¬ 
sider for instance the case g = 1/5. Writing sioj) = 
F{u/T b), where F is a universal function, let 
us compare the contributions to the imaginary part of 
F coming from other intermediate states. (The contri¬ 
butions for the real part have similar properties). This 
can be done by using explicit formulas for form factors 
deduced from |l^ . We display in figure 1 the next domi¬ 
nant contributions, i.e. the one particle form factor of the 
breather n = 3 and the two particle soliton/anti-soliton 
form factor. It is clear from the figure that these contri¬ 
butions, which can barely be seen on the graph, are negli¬ 
gible when compared to the one breather contribution for 
most of the range. In the deep IR, when w/Tg > 10, then 
the soliton/anti-soliton form factor starts to contribute 
(indeed it is necessary to get the correct IR exponent) but 
the function F is very small in that limit and this is cer¬ 
tainly not the dominant region. The perturbative expres¬ 
sion is also given on the graph to give an idea of the differ¬ 
ences. The g dependence of the correlator is then entirely 
included in the reflection matrix of the breather (p^) 


and is very easy to extract from the previous formulae. 



FIG. 1. Accuracy comparison for various form factor con¬ 
tributions for g = 1/5. 
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